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Abstract 

The mathematically correct computation of the spin glasses free energy in the infinite range limit crowns 
25 years of mathematic efforts in solving this model. The exact solution of the model was found many 
years ago by using a heuristic approach; the results coming from the heuristic approach were crucial 
in deriving the mathematical results. The mathematical tools used in the rigorous approach are quite 
different from those of the heuristic approach. In this note we will review the heuristic approach to spin 
glasses in the light of the rigorous results; we will also discuss some conjectures that may be useful to 
derive the solution of the model in an alternative way. 
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1 Spin glasses 

Spin glasses are interesting for many physical reasons; moreover they are the prototype of a complex system, 
i.e. a system with many phases. Here we will consider only spin glasses with an infinite range interaction, 
where the mean field approximation is correct (more of less by definition). In this case [TJ [3l |H [5] (i.e. the 
Sherrington Kirkpatrick model), the Hamiltonian is given by 



HM = 2 



r ^2 J^kViVk j (1) 



i,k=l,N 

where the N variables a are Ising spins, i.e. = ±1. The J are random variables: they are Gaussian 
distributed with zero average: 

E[J hk } = 0, E[jf >k ]=N- 1 , (2) 

where E[-) denotes the expectation value respect to the random couplings J. 

For each instance of the system (i.e. a choice of the couplings J) we can compute the ground state energy 
density that is defined as 

N _mm a Hj[a} 
ej ^ • (3) 
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The normalization factors have been chosen in such a way that the energy density is a number of order 1, 
when N — > oo. In the same way for each choice of the variables J we can define the free energy density 
fj (0) and the partition function Zf ((3) as 

f?W = - logi % m , zm = E eM-PHM) ■ (4) 

cr 

When N — » oo, it can be proved that the fluctuations in the energy density and in the free energy density 
go to zero with N: 

E[(em-(E[a»]) 2 ^0, E[(f? ) 2 ] - (E\f?])* - . (5) 

Therefore, with probability 1 when N goes to infinity, all systems have the same free energy density. 
We want to compute 

/(/?) = lim E[ff{p)] , g = lim E[ef] . (6) 

In the limit N — » oo the probability distribution of fj((3) becomes a delta function and the average value 
becomes equal to the most likely value. 

The limit N — > oo and /3 — > oo can be exchanged^ and therefore we have that 

g = lim /(/?) . (7) 

The mean values are not the only interesting quantities. We would like to control the whole probability 
distribution of the free energy (/) at fixed and large N (i.e. Pn(/) El)- We can consider various problems: 

• The behaviour of the typical fluctuations around the average value. For example it has been conjectured 
that in the region where the probability is concentrated we have [6j [7] : 

F w (/)«iV- 5/6 p(^ 6/5 (/-/)) , (8) 

where p is an unknown function and the difference between the r.h.s and the l.h.s. is negligible for 
large N (unfortunately we do not know how to compute heuristically the function p) . 

• Large deviations. For values of the free energy that is smaller than the most likely value (i.e. / < /), 
the probability is exponential small in N [BJ[I5]; the following limit exists (for not too large /): 

N~>oo 1\ 

• Very large deviations. For values of the free energy that is bigger than the most likely value (i.e. 
/ > /), it has been conjectured [HIS] that the probability decrease faster than any exponential of N 
and that it is exponential small in N 2 [7] . The following limit should exist (for not too large /) and it 
should be non-trivial: 

C{ f) = lim Z^0W1 . (10 ) 

Similar questions are well posed also for the ground state energy e. 

As we shall see below, there are rigorous explicit expressions for / [3] and for the large deviations function 
L(f) 10J . On the contrary for the very large deviations function C(f), there are only heuristic evaluations 
[7]. Nothing analytic is know (only numerical evaluations) on the function p that enters in the scaling form 
of the probability Pn(J) (eq. ([8])). The form in eq. J5|) has been conjectured in order to match with the 
known behaviour of the large deviation function, i.e. 

i(/)«(/~/) 6/5 (11) 



1 As suggested by a referee, the exchange of limits is permitted by the elementary inequality \E[fj (/?)] — E[e]y]| < ln(2)//3. 
2 From time to time, when this does not lead to ambiguities we will not indicate in an explicit way the /3 dependence of the 
various quantity: the correct, but heavier, notation would be Pjy(/, /3). 
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for / near (and smaller than) / and with the supposed behavior of the very large deviation function, i.e. 

^(/)«(/-/) 12/5 • (12) 

The reader should notice that the same kind of questions have a known answer for the ground state 
energy in the linear case where the variable a are continuous variables: 



N 

i=l,N 



Indeed in this case the ground state coincides with the largest eigenvalue of a Gaussian random matrix where 
nearly everything is known [111 112) . Here we are interested in the more difficult case, 

Gi = ±1 , (14) 

where ej is not a function of only the eigenvalues of J. 

2 Some heuristic considerations 

Let us start with a simple theorem |13j that gives interesting information on what we should not assume. 
For a given system we introduce the connected correlation functions at temperature /3 

Ci^k — (c.jCTfc} Qibbs — {o~i) Gibbs {o~k) Gibbs , (15) 

where the bracket denotes the average with respect to the Gibbs distribution^: 



(A[a]) Gibbs = lim £ P[a]A[a] , P[a] = ex P(-^H) . 



(16) 



It would be natural to assume that the connected correlations go to zero when N — > oo, more precisely 
to assume that 

V C 2 

lim ^ l ' k = o . (17) 

AT^oo N 2 

If the previous assumption would be valid, we could derive a simple formula for /(/3) that implies that 

^2 



e= lim /(/3) = -W-. (18) 

/3^oo V 7^ 

The results is reasonable, but in definite disagreement with numerical estimates for N up to O(10 3 ). However 
the situation is much worse: the same argument would imply that the entropy s(f3) satisfies the relation 

lim s(/3) = < . (19) 

P^oo (.271") 

The entropy s((3) is non-negative by definition. We can thus conclude by absurdum that 

lim ± . (20) 

N— too iv 

In order to understand the importance of the result (|20[) . it is convenient to recall some know facts on 
clustering states [1411151 [TB] . Let us consider an infinite translational invariant (non-random) system. In this 



3 In order to lighten the notation I have not indicated the dependence of both the Gibbs expectation values ((•)) and of the 
P[u] on /3 and J. 

4 The previous assumption should be formulated more carefully in order to take into account the invariance of the Gibbs 
measure when we change sign to all the spin simultaneously; I will not discuss here this technical detail. 
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case one can introduce the general concept of an equilibrium state, i.e. a state that satisfies local equilibrium 
conditions (e.g. the DLR equations 16J). If the Gibbs probability distribution has a limit when the volume 
becomes infinite, the Gibbs probability distribution is an example of equilibrium state, but it may be not 
the unique equilibrium state. 

Now the set of equilibrium states is convex; by the Alaoglu theorem any state can be written as the linear 
sum of the extremal (pure) states. Therefore we can write 

{■)Gtbbs = y^w n (') a , y] w a = i , (21) 

a a 

where the sum runs over the pure states (for simplicity I am assuming that the relevant set of pure states is 
discrete). A well-known theorem states that in a pure state connected correlation functions vanish at infinity 
(at large distances), i.e. pure states are clustering states, and vice versa. Intensive quantities (i.e. average 
on the whole system of local quantities) do not fluctuate in pure states. 

The ferromagnetic Ising at zero field at low temperature provides a well-known case: there are two phases 
(+ and -): 

(oi)± = ±m, (<Ji<Jk)± ~ m 2 , w + =w_ = i, (^Gibbs = to+(-)+ + to_(-)_ , (22) 
where, for i and k far away, 

{&i<J k) Gibbs ~ (o») Gibbs (o"fc) Gibbs = m 2 ; ((7j(7 fc } + — {(T l ) + ((7 k ) + W . (23) 

The situation is more complex for disordered systems, where the Gibbs measure may not have a limit 
when the volume goes to infinity. Although the appropriate and sophisticated tools have been forged to 
deal with this situation, it may be convenient to consider what happens in large, but finite systems, without 
putting the <5's and e's needed to make the statements sharp0. 

Let us consider a system where the number of degrees of freedom (TV) is large, but finite. We suppose 
[T71 [TS] that for large N we can write 



/Gibbs 



The states are approximately clustering if 



Hi,k^i,k) 

N 2 



o(l) . (25) 



The previous theorem (see eq. (|2U)) ) tells us that for spin glasses more than one w should be different 
from zero and therefore the decomposition in eq. (|2"4"|) is non-trivial. In order to describe better the Gibbs 
decomposition is convenient to introduce additional concepts. 

The first interesting quantity is the probability distribution of the w's (i.e. V({w})); however the weights 
w's of the states do not convey the whole information. We would like to know how the states differ one from 
the others: at this end we can introduce the distance and the overlap between two states: 

_ £i=i,jv (MO)a - MOM 2 _ Ea^jvHOM "^)^ r9R , 

«ct,7 — j 9a, 7 — ■ (,2b) 

Obviously the distance d and the overlap q are related. We have 

d an = 2(q EA - q a ,~i) , (27) 



5 As a physicist I would also remark that all the experimental systems are finite, and a statement on the behaviour of large 
finite systems can be experimentally tested much better than a statement on the behaviour of an infinite systems. 
6 Morc precisely more than two uj's. 
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where qEA = qa,a — 97,7 (we are assuming that q a . a is independent form a). Usually in literature one 
studies the properties of the overlap, the distances have been introduced here because their definition is 
more natural. 

Now for a given system J, we introduce the descriptor T>j, that is given by the following set 

} , (28) 

where the index a runs over a numerable sefl 

The descriptor Vj does depend on J. The functional V(V) is the probability distribution of the de- 
scriptors. One can construct heuristic arguments pQ to compute the probability distribution V(T>), the free 
energy / and the large deviation function L(f). We will describe these results later on in section (4). We 
now proceed into the description of some rigorous results. 



3 Rigorous results 

In this section we will concentrate our attention on some theorems. They are obtained using two different 
approaches, the first is based on probabilistic techniques, the second heavily relies on Guerra interpolation 
technique. 

Without any reference to the physical meaning of the descriptor and to the Gibbs phase decomposition 
we can introduce the descriptor T> that is defined as a normalized vector and a matrix, satisfying certain 
inequalities [19]. Neglecting technical details, a descriptor (T>) is given by the set of w a and of <? Q , 7 , where 
a belong to a numerable set and J2 a w a = 1. In this context one would naturally associate to a descriptor 
a decomposition of the Gibbs measure into states, however for the sake of proving theorems one can define 
a descriptor in an abstract way. Moreover we can also introduce a probability distribution on the space of 
the descriptors (V(T>)). The following theorems hold. 

• One can define a free energy functional T\P\ ('P(X') being the probability of descriptors) such [TH] that 
the true free energy is given by 

f{(3) = maxJ"[?] . (29) 

The explicit form of the free energy functional T\P\ is rather complex [19] and it will be not described 
here. 

• Using arguments, which do not refer to descriptors and are based on Guerra interpolation, it is possible 
to compute in an explicit way the value of / and consequently (from the previous theorem) max-p T\P\ . 
By inspection one verifies that the maximum is reached for the probability distribution predicted by 
the heuristic approach (V*): also the value of the free energy /(/3) coincides with the one computed 
by the heuristic approach [201 HI]. Therefore we know a probability distribution that maximize T\P\. 

• Under some assumptions, i.e. if stochastic stability and ultrametricity (see the definitions later) hold 
|22) . the function 

P(q) = *^2w a w~ ( 5(q - q ar( ) (30) 

determines the probability distribution ViJ)). Therefore in the stochastic stable ultrametric case the 
probability distribution is essentially unique, apart from a redefinition of the overlap. 

The bad new is that we do not know if the true V(T>) has some relation with the V*(T>) that maximizes 
the free energy; for the moment the probability distribution V(T>) does not have a definite meaning. It also 
not clear if the maximum of the free energy functional is unique, also if we consider only those probability 
that are stochastically stable H. 

7 The numbers of states is obviously finite for finite N, however in order to consider the limit N — » oo it is useful to consider 
a numerable set of states. 

8 The reader can find the definition of stochastic stability and a discussion on its consequences in |23l 1241 l25l I26| . 
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Figure 1: A taxonomic tree: the leaves denote the states: they are a the bottom. The distance among the 
states corresponds to the minimum height we have to cross for going from a state to an other state. This 
definition of the distance is ultrametric. 



It may be interesting to note the maximum is reached for the probability an ultrametric V*: i.e. the 
probability matrix is concentrated on distances that satisfy the ultrametric inequalities: 



Ultrametricity implies that the probability distribution of the distance among three random configura- 
tions is supported, in the limit of very large systems, only on equilateral and isosceles triangles with no 
contributions coming from scalene triangles. In an ultrametric space if two spheres have one point in com- 
mon, their union coincides with the sphere having maximum radius: in other words a random walk with 
steps of length I would never arrive at a distance larger than I from the origin. 

Ultrametricity is a very striking property: it implies that the equilibrium configurations of a large system 
can be classified in a taxonomic (hierarchical) way (as animals in different taxa). When ultrametricity holds, 
the space of the descriptors reduces to the space of infinite weighted trees (each leave has a weight). Ruelle 
gives a precise mathematical definition of a probability distribution over these trees [271 [28, 29 (in his paper 
one can also find interesting comments on the pruning of infinite trees.) 

3.1 Some details on the solution 

In this section we present the explicit form for the free energy f{(3). We firstly introduce an non-decreasing 
function q{x) that is defined in the interval [0 — I]. In the replica approach the inverse function q[x] is related 
to the function P(q) (defined in eq. ([30)) ) by the equation: 



Jo 

Let us consider a function g(x,h), that is defined in the strip < x < 1. The function g satisfies the 
following boundary condition at x = 1: 



Moreover (if we assume for simplicity that q(x) is differentiable), the function g(x,h) must be the unique 
solution of the following antiparabolic equation 



< max(d aj/ 3, d/3 i7 ) V/3 . 



(31) 




(32) 



<?(l,/i) = log(cosh(/%)) . 



(33) 




(34) 



G 



We now introduce a functional F[q]; 

F[q}=(3 f dx{\~q{x) 2 ) -5(0,0) . (35) 
Jo 

Talagrand theorem tell us that the free energy F is given by 

F = max F[q] . (36) 

q(x) 

Moreover, using Ruelle construction |27j . we can associate to each function q{x) a probability distribution 
over the descriptors, i.e. V{T>). For the probability distribution over the descriptors V(T>) that corresponds 
to a function q(x), the functional T\P\ (introduced in [T!|]) gives the same result of eq. ([33)1 . 

These results extend to large deviations. Indeed a theorem of Talagrand states that for 1 > n > we 
can define in the large N limit a large deviation function f(n): 

dP N (f) exp(-Nnpf) « exp(-Npf(n)) . (37) 

The functions f(n) and L(f) are simple related by a Legendre transform at usual. Talagrand has proved 
that the explicit value of f(n) is given by 



f(ri) = max 

q(x) 



[3 I dx(l-q(x) 2 ) -g(n,0) 



n < 1 q(x) for n < x < 1 . (38) 



It is evident that for n — we have that /(0) = /, as it should be. 

It is possible to perform analytic and numerical evaluations of the function f(n) following [ElOU], where 
the properties of the function f(n) are discussed in details. 

When n < 0, there are no theorems. There is a conjecture [5]: 

log(fdP N (f)exp(-Nn/3f)) 
ton BU U > P ^ = / = / . 39 

N^oc Nnp 

In other words the probability Pn(J) should be smaller than an exponential of N if / > /. 
What about very large deviations? In [7j it has been conjectured that 



dPj 



N 



(/) exp (-N 2 np{f - /)) « cxp (-N 2 0f LL (n)) . (40) 



Moreover some heuristic computations of the function jhhin) have been performed and one arrives the 
following formula: 

P N (f) « exp (-CiV 2 |/ - /| 12/5 + ...), (41) 
where the value of C is approximately evaluated [7] and it should be near to 2/3. 



4 The algebraic replica approach 

In the original heuristic approach the computations were done using the algebraic replica method [T] . Later 
it was found that the replica method was equivalent to a probabilistic approach [T] that inspired the rigorous 
results [501 dHl HI] • However it may worthwhile to try to understand the basis of the replica method in order 
to see if we can give it a rigorous basis. This is not easy because the algebraic replica approach is based on 
very strange mathematics. 

The mean field theory was solved by replica approach. The basic idea is quite simple and essentially it 
goes back to Nicola d'Oresme. Our goal is to compute 

/» = jU/3) = lim f N {J3) , f N (j3) = . (42) 
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However in our case it is much simpler to compute (for integer n) 



/ArK/?) - (3N7r ' (43) 



and at the end performing the limit 



foo(n,P)= lim f N (n,(3). (44) 



Although the computations are simple for integer n the construction can be done also for non-integer n and 
we have that 

Um /(n,/3)=/oo(/3) ■ (45) 

n— >0 

The function /^(n, /3) is interesting because it is the one relevant in the large deviations regime. 

A simple computations can be done for integer n, while we need the function /^(n, /3) also for non- 
integer n. While /jv(n., 0) is an analytic function n, at low temperatures (high (3) foc(n, (3) is not an analytic 
function of n so that we cannot obtain its value at n = by analytic continuation in n from integer n. 

In order to bypass this problem physicists were lead to an unusual mathematics: as we shall see below 
one introduces an n x n matrix Q, where eventually n is analytically continued to 0. 

An elementary computation (based on Gaussian integrals) tell us that when n is integer 

exp(-0nNf N (n,p)) = log (e [(Z£(/9)) n ]) = J dQexp(-0NnFf,[Q]) , (46) 

where Q are n x n symmetric matrices (that are zero on the diagonal) and the integral is done over all these 
matrices (the function F[Q] has a simple expression). 

It is evident that for positive integer n the integrals in equation (|46p can be evaluated by the saddle point 
method and one gets. 

/oo(n,/3) = min^IO]) = F [Q*] . (47) 
w 

Unfortunately, at high the function /oo(ra) that is defined for all n, is non-analytic in n. So in order to 
evaluate it in the interesting region and to find out the value of Q* one has to consider integral representation 
for non-integer n. In this way one is lead to consider the saddle point over 0x0 matrices. 

There are many ways in which a x matrix may be constructed. At the end of the day in the replica 
approach the matrix Q is parameterized in terms of a function q(x) defined on the interval 0—1, so that 
the space of x matrices becomes infinite dimensional. The game consists in computing everything by 
analytic continuation of the integral representation. 

There are many way equivalent to do the construction and they lead to the same result. We choose one 
that is particularly simple [30J, [31] . 

Let us take a prime number p and let us assume that there is an integer m such that p m — n. We look 
for saddle point of Fp[Q], where the matrix Q has the form: 

Q a ,h = q{\a-b\- x ) . (48) 

Here \a\ p is the p-adic norm: (the p-adic norm is defined as follows: if p k divides a and p k+1 does not divides 
a, \a\ p =p~ k ). 

We write Fp[Q] as function of p and m we make an analytic continuation in p and m of the result. 
Everything is computed by making analytic continuations from the region where the construction make 
sense. We will also write down the equations for the critical points in an explicit way for the function q(x). 
The final recipe consists in assuming that the free energy is given by the analytic continuation of the free 
energy that is evaluated at the analytic continuation of the saddle point. There are no serious mathematical 
justifications for this procedure. The heuristic justifications are: 

• At the present moment in the framework of the replica approach there is nothing else that you can try. 
You have two choices: or you accept this approach, or you change framework (or you invent something 
new). 
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• It works, it gives results that have sense and that have been confirmed by a rigorous mathematical 
analysis. 

• It is quite likely that the rigorous results would have not been derived without the heuristic analysis 
of the replica method. 

Now we can proceed to the limit n — > 0: at this point nothing is anymore an integer. Let us suppose 
that the prime number p satisfies the condition < p < 1 ; this is quite possible because p is not anymore a 
prime number, but it is the analytic continuation of a prime number. Now if we send m — > oo the quantity 
n goes to zero, as it should be (p is less than 1). 

As a final step we perform the limit p — > 1~: the range of \a — becomes the interval[ — 1] and 
the matrix Q depends on the function q(x) defined on the interval [0 — 1]. In this way the function F[Q] 
is promoted to a functional of q(x) and one finds the previous formulae for F[q]. The replica approach 
gives the exact result for the free energy and the functions q(x) is the same of the one of the mathematical 
approach. It is remarkable that integers are an ultrametric space with respect to the p-adic distance and 
this ultrametricity property implies the ultrametricity of the distances in the descriptor. 

Later one it was a surprise to discover was that all the formulas of the replica approach can be translated 
into probabilistic statements pQ. The matrix Q of the replica approach is a very compact to code the 
probability V(T>). It looks strange that a matrix codes a probability distribution, so that let us explain how 
it works. 

One can prove the following formula 



E [S(Q a ,b - «)] = E 



w a wpS(q ai - q) 



(49) 



where the E at the l.h.s. denotes the average over the indices a and b (with a ^ b) and E at the r.h.s. 
denotes the average over the probability distribution of the descriptors. In a similar way one can prove 

E [6{Q a , b - qi, 2 )S{Q b ,c - q2,s)S(Qo,a ~ q 3 ,i)] = (50) 

E " 



, w a w !3 w 1 5(q a j j - gi, 2 )<%/3, 7 - <72,3)% 7 

,a 93, lj 

Going on one can write an infinite set of relations that should uniquely determine the probability distribution 
over the descriptors. It is quite remarkable that a x matrix is used to code the probability distribution 
over a space of infinite dimensional matrices and that the ultrametricity of the integers with respect to the 
p-adic norm implies the ultrametricity of the matrix of distances among states. 

The reader should note that not all Q matrices codes for a probability: the probability must a positive 
function and there are many matrices Q where the l.h.s. of eq. (|49[) cannot be interpreted as a probability. 
It is also not clear if there are probabilities V(T>) that cannot be coded using the matrix Q: a negative 
result that is difficult to prove and there is no constructive procedure to get the matrix Q (if any) from the 
probability distribution of the descriptors. 

In conclusions the heuristic version of the descriptor approach (i.e. the cavity method) and the replica 
approach are equivalent [17] ■ The functions q(x) of both approaches are equal. The replica approach gives 
the exact result for the free energy. 

The reader may be perplexed by this procedure that was invented about thirty years ago. However it 
has a very strong heuristic value and both the heuristic probabilistic approach and the rigorous results have 
strongly beneficed by the fact that the answer to many question was known especially in this context that 
is rather unfamiliar. 



5 Two conjectures 

Can the original replica derivation made rigorous? The distance from conventional mathematics seems to be 
very far so that a direct assault may be not successful. However it is clear that the subject of our discussions 
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is the set of critical points of the function F[Q] for all n, in the same way as the analytic continuation of a 
function on the integers is related to the value that the function takes on all integers. 

Following [32] we would like to present some conjectures, whose proof would be a step forward a rigorous 
understanding of the replica approach. 

Let us consider a function Z(n, N) that can be simultaneously written as 

r+ca /• 

Z(n,N) = J dz(i N (z)z n = J dQexp(—NF(Q)), (51) 

where (a(z) is a positive function and the integral runs over symmetric n x n matrices; F(Q) is an analytic 
function of n, in the sense that it has a fixed form for all n 0. The existence of functions, which can be 
written in both ways, is non-trivial, however the partitions function of spin glasses and generalized spin 
glasses provide an explicit example. 

We know define a new function Z(n, N), i.e. 

Z(n,AO = ]Texp(-iVF(Q fe )) (52) 

fe 

and the sum goes over all the critical points of the function F, i.e. those satisfying the condition 

OF 



dQ 



= 0. (53) 

Qk 



. It not clear which of the two definitions is the most appropriate. 

Qk 



An alternative definition could be Z(n,N) — J^k Sk ex P( — N F(Qk)), where the factors Sk are given by 
s k =sign (det(0)) 

Our aim would by to approximate Z(n, N) with Z(n, N) also for non-integer n. At this end we need to 
have some a priori control on the analytic properties of both functions. In this context we put forward the 
following conjectures. 

• The function Z(n, N) satisfies for integer n the following integral representation. 

r+oc 

Z{n,N) = / dzyi N {z)z n . (54) 



The previous formula implies the existence of a natural (and under certain conditions unique) analytic 
continuation of Z(n, N). 

The function Z(n, N) is a good and uniform approximation to Z(n, N) for large N: we should have: 



lot 



(z{n,N) 



km >- = km kg (55) 

N^oo N AT-i-ao N 

The previous formula should not only for integer positive n, where it is trivial, but also for non-integer 
n and in particular for n in the interval [0 — 1]. In the best of possible words one should have: 

. lim MM**» 

JV^oo N JV—oo N 

These conjectures are interesting also because not many results are known in this direction (as far as I can 
tell). However in the contest of this paper they are very relevant because in some cases these conjectures 
have been used to prove the result of coming from the replica approach [32) . 



9 Examples of allowed functions are Tr(Q 4 ) or . the factor 1/n is needed in order to have a non zero limit at 

n = in interesting cases. 
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